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We construct a quantum extension of the (classical) three-coloring model introduced by Baxter [J. Math. 
Phys. 11, 784 (1970)] for which the ground state can be computed exactly along a continuous line of Rokhsar- 
Kivelson solvable points. The quantum model, which admits a local spin representation, displays at least three 
different phases; an antiferromagnetic phase, a line of quantum critical points, and a ferromagnetic phase. We 
argue that, in the ferromagnetic phase, the system cannot reach dynamically the quantum ground state when 
coupled to a bath through local interactions, and thus lingers in a state of quantum glassiness. 



I. INTRODUCTION 

It has long been known that strongly correlated electronic 
systems at zero temperature can be parametrically close to 
fixed points with macroscopically degenerate ground state 
manifolds. This is so for a two-dimensional degenerate elec- 
tron gas in the quantum Hall regime for which the largest 
energy scales are the cyclotron and Zeeman splittings, and 
for Mott insulators when the Coulomb energy in the form of 
an on-site repulsive Hubbard U dominates all other energy 
scales. 1,2 Similarly, the zero temperature phase diagram of a 
class of quantum dimer models and geometrically frustrated 
quantum systems are known to be influenced by fixed points 
with massively degenerate ground state manifolds. 3 Whereas 
most of the efforts have been directed toward exploring the 
nature of the collective excitations present in the different 
phases and at their boundaries, relatively little attention has 
been paid to nonequilibrium phenomena induced by the cou- 
pling to a bath. The goal of this paper is to illustrate how 
quantum glassiness can emerge in a two-dimensional homo- 
geneous quantum system whose local degrees of freedom are 
subjected to a strong constraint, when the system is cou- 
pled to a bath. In the process, we also show how to obtain 
a whole line of solvable Rokhsar-Kivelson (RK) points in a 
hard-constrained quantum system. 

More specifically, we are going to construct a quantum 
version of the (classical) three-coloring model introduced by 
Baxter in Ref. 4 which we shall call the quantum three- 
coloring dimer model. We find that: (i) For any value of a 
dimensionless coupling (3 J the ground state follows from tak- 
ing a superposition of states in a preferred basis of the Hilbert 
space. The dependence on (3 J of the normalization of the 
ground state defines the partition function of a generalization 
of the classical three-coloring model at the reduced inverse 
temperature f]Jm&^ (ii) The excitation spectrum above the 
ground state is identical to the spectrum of relaxation times of 
a classical rate equation that depends on a dimensionless pa- 
rameter a in addition to f}J. (iii) The ground state expectation 
value of any two local operators at unequal points in space and 
imaginary time can be interpreted as the classical correlation 
function of the classical three-coloring model endowed with 
dynamical rules, provided both operators are diagonal in the 



preferred basis. This interplay between a quantum model and 
a classical stochastic model generalizes a similar connection 
established by Henley for the square lattice quantum dimer 
model. 8 (iv) The quantum three-coloring dimer model is lo- 
cal in that it can be interpreted as a Hamiltonian for quantum 
spins whose interactions have a. finite range as long as the pa- 
rameter a remains nonnegative. In the spin representation, the 
zero-temperature phase diagram has at least three phases, an 
antiferromagnetic (AF) phase when — f3 J 3> 1, a line of quan- 
tum critical points in the neighborhood of infinite temperature 
and a ferromagnetic (F) phase when (3 J 3> 1. (v) We argue 
that, in the F phase, the relaxation time needed for a generic 
initial state to decay to the ground state through a local cou- 
pling to a bath diverges exponentially with the system size, a 
property that signals quantum glassiness^ 



II. DEFINITION OF THE QUANTUM THREE-COLORING 
MODEL 

The classical three-coloring model is a model of dimers laid 
down on the bonds of a lattice with coordination number 3. 
More precisely, the dimers come in three kinds (or colors), 
A, B, and C and occupy the bonds (ij) between nearest- 
neighboring sites i and j of the honeycomb lattice according 
to two rules: (1) all bonds are occupied (close-packed con- 
straint); and (2) all three dimers meeting at a site (vertex) 
are different (three-coloring constraint). The classical phase 
space is thus the set S of all configurations C obtained by col- 
oring the bonds of the honeycomb lattice in any of the three 
colors A, B, C with the condition that three bonds meeting at 
a site have different colors. Two possible configurations C AF 
and C F are depicted in Fig.^ 

To any configuration C E S there corresponds a unique con- 
figuration {cr s } of Ising degrees of freedom a i = ±, defined 
on the sites i of the honeycomb lattice, whereby a { = +(— ) 
if the three dimers meeting at i can be obtained from ABC 
through an even (odd) permutation when moving clockwise 
around i. The configuration C AF {C F ) in Fig. [2 thus corre- 
sponds to an AF (F) configuration of the Ising spins. Evi- 
dently, any global even permutation of the colors ABC leaves 
the Ising configuration unchanged. Hence there is a three-to- 
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FIG. 1 : (Color online) Two configurations, Caf m panel (a) and C F 
in panel (b), from the three-coloring model. 



one relation between S and the set of all Ising spin configu- 
rations on the honeycomb lattice such that the magnetization 
along the perimeter of any elementary hexagon is either as 
in Fig.[3a) or ±6 as in Fig.[2b)>i£ 

Another useful characterization of a configuration C G S is 
in terms of two-color paths, each of which is defined by draw- 
ing a line along all bonds colored in any given two (say B or 
C) of the three colors that initiate from a site (the seed) and 
by specifying the color pattern (say BCBC... or CBCB...) 
along the line. Upon imposing periodic boundary conditions 
one can interpret a configuration C G S as a set of two-color 
closed paths (where the same two colors have been chosen 
for all the paths in the lattice) that are nonintersecting and 
close-packed. We will refer to such objects making up a con- 
figuration C G S as loops. For example, configuration C AF 
(C F ) in Fig.[Qcan be viewed as being made solely of loops of 
maximal (minimal) nonnegative curvature. 

We will denote by £ a decorated loop, i.e., an object made 
of all 2L e independent colored bonds emanating from the L e 
sites through which a loop passes, and by the set of its colored 
dangling bonds, as defined hereafter (see also Fig.|2j. For any 
£, the L e bonds that alternate between the colors B and C, 
say, form the backbone of I, while the remaining L t bonds of 
color A that stick out of the backbone connect the loop to its 
dangling bonds, i.e., the two bonds that depart from the free 
end of each A colored bond. The definition of a decorated 
loop therefore encodes the information about the colors of all 
the dangling bonds of the loop. The reasons for such a choice 
will become evident as we proceed with the formulation of 
our quantum model [see Ea.fTfl and thereafter]. The set of all 
possible decorated loops is C. Configuration C AF in Fig.^is 
made exclusively of objects I with six connecting bonds (i.e., 
the bonds joining the loop backbone to the dangling bonds) of 
color A, B, and C, respectively. Configuration C F in Fig. ^ 
is made exclusively of objects £ with L connecting bonds of 
color A, B, and C, respectively, where L is the linear size of 
the system. 

By specifying the colors and the locations of both the con- 
necting and the dangling bonds, one identifies two decorated 
loops £ and £ related by the color-swap operation B <-» C 
along the L e bonds making up the loop backbone. This color 
swap relating two different configurations C and C in S is il- 
lustrated in Fig. |2 Given a decorated loop £ and a configura- 



tion C, it is always possible to decide if £ belongs to C, i.e., 
either I € C or £ £ C. By construction: (1) If £ G C then 
£ £ C. (2) If I G C then there exists a unique configuration 
C (to which I belongs) obtained by the color swap operation 
along the backbone of £ in C. 




FIG. 2: (Color online) Two examples of decorated loops of the types 
BCBC... and CBCB... in the three-coloring model are shown. 
They are related by the color-swap defined in the text. Notice the in- 
formation encoded in both loops, (1) the backbone structure (circled 
lattice sites) that is specified by all the sites visited by the sequence of 
bonds of color B or C emanating from a seed, (2) the two-color pat- 
tern (BCBC... and CBCB...) along each loop, and (3) the colors 
and locations of all the dangling bonds (thick lines). The connect- 
ing bonds that link the backbone of a decorated loop to its dangling 
bonds are represented by wavy lines for clarity. 

To define the quantum three-coloring dimer model we need 
a Hilbert space and a Hamiltonian. The Hilbert space is the 
span of all orthonormal states \C) labeled by the classical con- 
figurations C G S. The Hamiltonian, which we show below to 
be of the RK type, takes the form 



eec 



(la) 

Q e ■■= v e - T tJ , (lb) 
where the potential and kinetic energy operators are 

V e := e-^^Pf, = e"^/ 2 ^ c e \C)(C\, (lc) 

ces 



T - ■= 



ces 



c e \C){C\, 



(Id) 



respectively. The dimensionless coupling f3J is real valued 
and c t = 1 if £ G C, while it vanishes otherwise. The dimen- 
sionless energy 



(le) 



(ij) 



is the classical energy of the Ising spin configuration {aA- as- 
sociated to C in units of |J|. Finally, the dimensionless energy 



(If) 



is uniquely defined for any C G S that contains £. Notice 
that the definition of a decorated loop was precisely chosen 
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so that the energy difference — E c depends only on the 
information contained in I, and not on further details of C and 
C. 

Whatever the choice made for the energy scales {w e > 0}, 
the wave function 

l*^) = E ex P(-' 3Ji? c/2)|C) (2a) 
ces 

is the ground state of H RK , as follows from the facts that: (1) 
Q\ is proportional to Q e , and (2) \^ RK ) is annihilated by Q 
(i I E C). Remarkably, the normalization of 
nothing but the partition function 



rk&rk) is 



Z:=J2e- 

ces 



PJE a 



(2b) 



of the generalized three-coloring model (see Refs. HE) with 
the configuration energy JTeb . i.e., up to a multiplicative factor 
of 3, it is the partition function of the nearest-neighbor Ising 
model on the honeycomb lattice with the local constraint that 
the magnetization on any elementary hexagon is or ±64£ 
The zero-temperature phase diagram of the quantum three- 
coloring model along the axis f3 J £ K contains the phase di- 
agram of the generalized classical three-coloring model along 
the reduced temperature axis J S K. 



PHASE DIAGRAM OF THE CLASSICAL 
THREE-COLORING MODEL 



III. 



The partition function i2b\ was evaluated by Baxter at in- 
finite temperature. 4 Through a height representation of the 
three-coloring model, it was shown in Ref. [H] that at infinite 
temperature the probability that two points of the honeycomb 
lattice separated by a distance r belong to the same loop de- 
cays with r algebraically and that this decay is captured by the 
811(3)! Wess-Zumino-Witten (WZW) theory. The finite tem- 
perature phase diagram was studied numerically in Ref. 5 with 
a transfer matrix approach, in Ref. |g with a cluster variational 
method (CVM), and in Ref. 7 using Monte Carlo simulations, 
the CVM, as well as perturbative renormalization-group (RG) 
techniques for infinitesimal values of the coupling (3 J. The 
classical configuration that minimizes the energy E c is the 
antiferromagnatic C AF when J < and the ferromagnetic 
C F when J > (see Fig. [Q. Thermal fluctuations are ex- 
pected to destroy both long-range ordered phases for a suf- 
ficiently large temperature. In view of the local constraint 
of the three-coloring model, the transition may be different 
in nature from the one of the unconstrained model. Numer- 
ical simulations are consistent with the termination of the F 
phase in a first order phase transition at the reduced temper- 
ature (/3 J) F i On the other hand, the 811(3)! WZW critical 
point at infinite temperature seems to be the end point of a 
line of critical points when < 0J < ((3J) Qr . Whether 



{PJ) F or there exists intermediary phases aii in 



the range (/3J) CT < f3J < {/3J) F is not known at present. 
The phase diagram is not symmetric under J — > — J as the 



local constraint favors the AF over the F state from an en- 
tropic point of view. The CVM predicts a critical end point 
{/3J)ap < to the AF phase followed by a paramagnetic 
phase, possibly of the short-range resonating-valence-bond 
(RVB) type, that terminates in the critical point at infinite tem- 
perature. The conjectured phase diagram of the partition func- 
tion ( I2bl is sketched in Fig. [3] 



AF 



SU(3) 



FIG. 3: (Color online) Conjectured phase diagram of the classical 
partition function J2bt along the reduced temperature axis j3J . The 
right end point of the AF phase is a hexagon. The left and right end 
points of the critical phase are a circle and a square, respectively. The 
circle corresponds to /3 J = and realizes the SU(3)! WZW critical 
theory. 



IV. EXCITATION SPECTRUM, LOCALITY, AND 
EMERGING TOPOLOGICAL QUANTUM NUMBERS 

Although the ground state of H RK does not depend on the 
choice of the couplings w g , the excitation spectrum does. We 
are now going to show that there exists a one-to-one corre- 
spondence between the excitation energies of H RK and the 
relaxation modes of the three-coloring model endowed with 
the dynamics borrowed from a uniquely defined transition ma- 
trix (W cc ,). We begin by defining the matrix elements 



W, 



cc 



-P J ( E c- E «)/*{C\H RK \C'). (3a) 



Although (Wcc) i s not symmetric in general, it shares the 
same set of eigenvalues as H RK . Furthermore it satisfies by 
construction the detailed balance condition 



W cc ,e-P JE o = W c , c e-? JE c 



as well as the condition 



0= E W cc- 

ces 



(3b) 



(3c) 



Equation d3cl allows one to interpret the first order linear dif- 
ferential equation in imaginary time r 



Pc( T ) = Y] w ccPc( T ) 



ces 



(4) 



as a master equation for the probabilities {Pci T )} whose sta- 
tionary solution is 



pf = Z- 1 eM~(3JE c ) 



(5) 



by Eq. Obt . Furthermore, the properly normalized ground 
state expectation value G ~ ~ (r) of the product B(t) x A(0) at 
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unequal imaginary times of any two operators A and B diago- 
nal in the basis \C) can be interpreted as a classical dynamical 
correlation function at unequal times: 

d §A^ = Y,Y, B c>Pc'\c(t)A c p { ° } ^G ba {t). (6) 



The conditional probability p c , , c (r) is the solution to the mas- 
ter equation (0} with the initial condition p c ,\ c {t = 0) = S c , c 

while B c , and A^are the matrix elements of B and A, respec- 
tively. Equation l(6j, relating classical stochastics of the three- 
coloring model to the quantum dynamics in imaginary time 
of the quantum three-coloring dimer model, is an example of 
a generic property of quantum dimer models at the so-called 
RK point as noted by Henley. 8 

We can now motivate our choice for the couplings w e , 

w i = w ofi ex P i~ a ( L i " 6)] ■ (7) 

Here, the energy scale w Q is of order 1. The number < 
f e <l can always be chosen such that the classical stochastics 
in time encoded by Eq. is the Glauber or Metropolis dy- 
namics. 14 We do not believe this choice to be of consequence 
for spectral properties beyond the mean level spacing. The 
choice for the value of the dimensionless parameter a can af- 
fect spectral properties in a crucial way, however. Any a > 
implements (defines), in the thermodynamic limit, the condi- 
tion of locality. A corollary to a > is the emergence of 
well-defined topological sectors in the thermodynamic limit. 

The quantum three-coloring Hamiltonian Q where the 
coupling constants w e are defined in Eq. @ with a > is 
local in very much the same way as the effective quantum 
spin- 1/2 model obtained from the expansion of the Hubbard 
model at half-filling, and in the limit of an on-site repulsion U 
much larger than the band width t. The potential energy ( fTcl 
associated with a decorated loop I with perimeter L e generates 
2, 3, L £ -body Ising spin interactions for spins that can be as 
far as L e lattice spacings apart very much in the same way as 
ring exchange interactions along a closed path of length L g are 
induced to L £ -th order in t/U for the Hubbard model at half- 
filling. In the case of both the quantum three-coloring and the 
effective Heisenberg model, these L^-body interactions are as- 
sociated to characteristic energy scales that are exponentially 
suppressed with L e ; namely, w e and (t/U) L '- x t, respectively. 
The same argument applies to the kinetic energy dldt . In fact, 
we shall show shortly that here the counterpart to the Hubbard 
model is the constrained transverse field Ising on the honey- 
comb lattice. 

The concept of topological sectors is rooted in the follow- 
ing conservation laws characterizing any C G S. 4 Choose the 
brick wall representation of the three-coloring model as de- 
picted in Fig. [4] Draw a line between any two consecutive 
rows and count the number of intersections N A B c with ver- 
tical bonds of color A, B, and C, respectively. This triplet of 
numbers is independent of the choice of the two consecutive 
rows and characterizes globally each configuration C G S. 
Furthermore, for any I G C and for any C G S such that 
£ G C, N A B C = N ABC whereby N A B C (N A B C ) char- 
acterizes globally C (C), as long as t is nonwinding. The num- 
bers B c are elevated to the status of conserved quantum 
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FIG. 4: (Color online) Two configurations C AF in panel (a) and 
C F in panel (b) in the brick wall representation. In panel (a) iV^ = 
N B — N c . In panel (b) N A is maximum with N B = N c — for 
C F . 



numbers under the quantum dynamics defined by H RK , pro- 
vided w e = for all winding loops. In particular, these topo- 
logical conservation laws emerge in the thermodynamic limit 
when a > 0. It is possible that < j3J < (f3J) cr realizes a 
line of quantum Lifshitz critical points when a > 0, 15 with a 
dynamical exponent z = 2; this is in principle testable numer- 
ically. 



V. QUANTUM THREE-COLORING MODEL AND THE 
CONSTRAINED QUANTUM ISING MODEL IN A 
TRANSVERSE FIELD 

The RK three-coloring Hamiltonian fiat is a special case of 
the Hamiltonian 



eec 



(8) 



where the dimensionless couplings v g and t g are real valued. 
The choice of the couplings v e — ~e e J/(4:W e ) and t t = 1 
implies that: (1) the matrix element 

-(J/A)J2e e (C\PAC) 

eec 

is nothing but the classical Ising energy 

of the nearest-neighbor Ising model on the honeycomb lattice 
subjected to the local constraint on any elementary hexagon 10 



( 27T £ of/3 j = 1 

\ i£hex / 



(9) 



while (2) each term w l T l , is induced order by order in de- 
generate perturbation theory by the small transverse field 
r a f- Here, it is the local constraint that is responsible for 
the macroscopic quasidegeneracy of the classical limit V — > 
as opposed to the geometrical frustration in the Ising models 
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treated in Ref. This happens as long as the energy split- 
ting of order J characteristic of the classical configurations 
that obey the three-coloring constraint is negligible compared 
to the very large characteristic energy U paid by any classical 
configuration that violates this constraint. The energy scale 
associated to the three-coloring constraint plays a role simi- 
lar to that of a large Hubbard on-site repulsion for strongly 
correlated electrons or of the cyclotron energy in the regime 
of the fractional Hall effect. The constrained quantum Ising 
Hamiltonian in a small transverse field, 

^ing : = E (-V/fytPt - w eT itl ) , (10) 
lec 

captures the Z 2 quantum dynamics induced by time-reversal 
symmetry breaking in a triangular array of Josephson junc- 
tions when the superconducting order parameter breaks time- 
reversal invarianceiiii The quantum three-coloring Hamilto- 
nians H 3 _ c and H RK are also relatives to a model for quan- 
tum spin-1 degrees of freedom defined on the bonds of the 
honeycomb lattice introduced by Wen in Ref. 

We can now refine the notion of locality {a > 0) of H RK . 
The quantum dynamics in real time of H RK and of the con- 
strained quantum Ising model H lsj originate from the very 
same quantum tunneling processes between classical configu- 
rations. Their local nature is explicit in H lsing and carries over 
to H RK . 

VI. LOCAL COUPLING TO A THERMAL BATH, 
ARRESTED QUANTUM DYNAMICS OR QUANTUM 
GLASSINESS 

So far we have considered the quantum three-coloring 
model (JSJl or, equivalently, the constrained Ising model in a 
transverse field JlOb in isolation. Energy is then a conserved 
quantum number and the real-time evolution of any state in 
the Hilbert space is a linear superposition of harmonics with 
frequencies ojj = Ej/K where the index j labels the exact 
eigenvalues of Hamiltonians l|8} or dlOt . 

In the spirit of Caldeira and Legget, 19 imagine coupling 
H lsing (H RK ) very weakly (g lX < J, w ) to a bath of har- 
monic oscillators {a - A , a i X } through the local interactions^ 

S^of (oi.A, + h - c ) • (U) 
ix i 

The effect of this phonon bath is to endow the energy eigen- 
values of the strongly correlated system in isolation with finite 
lifetimes. These lifetimes can be calculated perturbatively. To 
put it differently, the presence of a phonon bath implies that 
we cannot consider a pure state, as was the case in deriving 
Eq. l|6}. The introduction of a density matrix is now required 
to study the thermodynamic behavior of H Ising (H RK ) and, in 

this sense, the evolution in real time of H lsing (H RK ) is better 
understood in terms of transition or relaxation rates between 



states of H IsiD „ (H RK ) by absorption or emission of phonons. 
At zero temperature, any initial state of H Ising (H RK ) which is 

different from the ground state of H uing (H RK ) can now decay 
to the ground state through spontaneous emission of phonons 
(aj \ ), i.e., by transferring energy to the thermal bath of har- 
monic oscillators. The goal of this section is to estimate the re- 
laxation rate of an initial state of H lsing (H RK ) into the ground 

state of H lsing (H RK ) due to the coupling to the bath. It should 
be emphasized that these relaxation rates are completely un- 
related to the imaginary-time dynamics of H RK encoded by 
Eq. ©. 

Deep in the F phase (J 3> w ) all the loops are straight and 
winding around the system (see configuration Cf in Fig.^. 
The relaxation time for an initial state to decay into the F 
ground state is bounded from below by a time of the order of 
the time needed by any state l^^gc ) = ^» f) to decay 
into the F ground state F ) through spontaneous emission of 
phonons. The matrix element for such a quantum tunneling 
process is obtained from the L-th order perturbation expan- 
sion of the coupling ( II It . and it is therefore suppressed by 
an exponential factor of order (g/U) L , where U is now the 
(large) characteristic energy scale for the three-coloring con- 
straint and g is the typical coupling to the bath [notice that 
this process requires 0(L) phonons]. Hence the relaxation 
time for an arbitrary initial state to decay through spontaneous 
emission of phonons into the F ground state is bounded from 
below by some positive power of exp(L). 9 The same rational 
applies to H RK . 

Exponential dependence of the relaxation times on the sys- 
tem size is a signature of the breakdown of thermodynamic 
equilibration (arrested quantum dynamics). For both H RK 
and -^1^, ^ occurs in the zero temperature phase diagram of 
Fig. [3] as soon as the ground and first excited states belong to 
different topological sectors. This is a direct consequence of 
the fact that changing the topological sector of a given state 
&)iLi requires flipping all the values of the a\ (via the 
corresponding a\ operators) along winding loops that have a 
length of order L. In the systems considered here, it so hap- 
pens that the change in topological sector of the ground state 
occurs only when the system orders ferromagnetically, i.e., for 
sufficiently large (3J . In fact, for all (3 J < {/3J) C1 . the ground 
and first excited states belong to the topological sector of the 
AF state, i.e., the typical loops that connect the ground to the 
first excited states have a perimeter of order 1, and no arrested 
quantum dynamics is to be expected. For sufficiently large 
(3 J the ground state belongs to the topological sector of the F 
state while generic excited states belong to different topolog- 
ical sectors and arrested quantum dynamics rules. 

An intriguing question is what is the fate of this arrested 
quantum dynamics when we forbid tunneling processes be- 
tween different topological sectors before taking the thermo- 
dynamic limit, and work exclusively in the topological sec- 
tor of the AF state. We believe that the mechanism for ar- 
rested quantum dynamics outlined above survives because the 
ground state for sufficiently large (3 J should have a finite over- 
lap with a polycrystalline F state, i.e., a state made of macro- 
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scopically large ferromagnetic domains in the Ising represen- 
tation or one made almost exclusively of very large loops in 
the loop representation (see Fig. 9 of Ref. 0). Indeed, this 
polycrystalline F state can be obtained from a generic excited 
state of the AF state (characterized exclusively by small loops) 
only by applying af operators along very large (though not 
system-spanning) loops, i.e., after the spontaneous emission 
of a number of phonons of order L. 

Dynamical arrest is one of the characteristics of glassy sys- 
tems. Classical glass formers, for example, fall out of equi- 
librium when the temperature is lowered and the relaxation 
times become of the order of the experimentally accessible 
time scales. In classical systems the fast decrease of the re- 
laxational rates with decreasing temperature is associated to 
the increase in the times needed to surpass energy barriers by 
thermal activation. In quantum systems, tunneling under high 
(but narrow) activation barriers remains an open channel for 
equilibration down to zero temperature. Such quantum tun- 
neling processes are the working principle behind the notion 
of quantum annealing, 20 which was demonstrated experimen- 
tally in the disordered Ising magnet LiHoo. 56 Yo.44F4^i How- 
ever, quantum annealing is rendered inoperative if the tunnel 
barriers are wide, and classical thermal activation may be a 
more efficient thermalization process even at rather low tem- 
peratures. In this case, quantum mechanics does not prevent 
systems from remaining in a state of glassiness at low temper- 
atures, very much as their classical counterparts can do in the 
presence of strong frustration induced by geometry or disor- 
der. 

What makes such a state of quantum glassiness remarkable 
is that it can occur in a system devoid of geometric frustration 
or disorder, solely as a result of strong constraints. At zero 
temperature, as thermal activation over energy barriers is shut 
down, the only mechanism for relaxation is quantum tunnel- 
ing. In the quantum three-coloring model, the effective tun- 
neling barrier widths are related to the sizes of the loops that 
are color-flipped. As the system tries to decay into polycrys- 
talline F states with larger and larger grain sizes, tunneling of 
larger and larger objects is needed, and hence the relaxation 
rates grow exponentially with the size of the loops as long 
as the quantum dynamics is of a local nature. The extreme 
case occurs when a generic out-of-equilibrium state attempts 
to decay into the F state, a process that requires the exchange 
of the colors along loops with a perimeter of the order of the 
linear size of the system. 

We have focused on dynamical arrest as a signature of 
glassy behavior. One may wonder whether a definition of 
a glass requires a thermodynamical transition into a glassy 
phase. The answer is negative, as many examples have been 
constructed of classical systems with trivial thermodynamics 
that display relaxation rates characteristic of glasses of the 



strong and fragile types. 22 Reference 9 provides one example 
of a quantum system with no thermodynamic transition but 
characterized by relaxation rates of a strong glass type; this 
system defies quantum annealing because the tunneling bar- 
riers become increasingly wide as the system effective tem- 
perature is lowered, and hence it remains in a state of quan- 
tum glassiness. The quantum three-coloring model is a second 
example thereof. In both cases, quantum glassiness emerges 
from the local nature of the Hamiltonian combined with the 
presence of strong constraints on the quantum dynamics. 



VII. CONCLUSIONS 

In this paper we have constructed a specific quantum exten- 
sion H RK to the three-coloring model of Baxter. Our strategy 
is very general as it applies to any space S of classical config- 
urations C and to any chosen subset C G S x S that defines 
which pairs of classical configurations are related by quantum 
tunneling. The strategy is specific in that it demands choos- 
ing (1) a potential term V{C), VC G S and (2) a kinetic term 
f(C, C), V(C, C) G C. Our choice for H RK was simple in that 
the ground state is known for all values of two dimension- 
less couplings a > and f3J G K. The zero-temperature 
phase diagram of H RK as a function of (3 J G M is intri- 
cate. Long-ranged ordered phases in the local Ising variables 
at (3 J <C — 1 and f3J ^> +1 are separated by qualitatively dif- 
ferent paramagnetic phases of which some are critical with re- 
spect to equal-time correlation functions. The excitation spec- 
trum of H RK depends crucially on whether a > (locality) 
or a = (nonlocality) and can be simulated using classical 
Monte Carlo methods. In the ferromagnetic phase, we have 
argued that the relaxation time for a generic state to decay into 
the ground state diverges exponentially with the system size, a 
trademark of quantum glassiness, when the coupling of H RK 
to a bath preserves the condition of locality. Thus the system 
prepared in a state out of equilibrium while in the ferromag- 
netic phase cannot decay into the exact ground state \^ RK ) 
by spontaneous emission of phonons from the bath and stalls 
into a mixed state however small the bath temperature is. 
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